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ABSTRACT 

Context. The thermal structure of a starless core is crucial for our understanding of the physics in these objects and hence for our 
understanding of star formation. Theory predicts a gas temperature drop in the inner ~5000 AU of these objects, but there has been 
no observational proof of this. 

Aims. We performed VLA observations of the NH3 (1,1) and (2,2) transitions towards the pre-stellar core L 1544 in order to measure 
the temperature gradient between the high density core nucleus and the surrounding core envelope. Our VLA observation for the 
first time provide measurements of gas temperature in a core with a resolution smaller than 1000 AU. We have also obtained high 
resolution Plateau de Bure observations of the 110 GHz 1 1 1 - loi para-NH^D line in order to further constrain the physical parameters 
of the high density nucleus. 

Methods. We combine our interferometric NH3 and NH2D observations with available single dish measurements in order to estimate 
the effects of flux loss from extended components upon our data. We have estimated the temperature gradient using a model of the 
source to fit our data in the «, v plane. As the NH3(1,1) line is extremely optically thick, this also involved fitting a gradient in the NH3 
abundance. In this way, we also measure the [NH2D]/[NH3] abundance ratio in the inner nucleus. 

Results. We find that indeed the temperature decreases toward the core nucleus from 12 K down to 5.5 K resulting in an increase 
of a factor of 50% in the estimated density of the core from the dust continuum if compared with the estimates done with constant 
temperature of 8.75 K. Current models of the thermal equilibrium can describe consistently the observed temperature and density in 
this object, simultaneously fitting our temperature profile and the continuum emission. We also found a remarkably high abundance of 
deuterated ammonia with respect to the ammonia abundance (50%+20%), which proves the persistence of nitrogen bearing molecules 
at very high densities (2 • 10^ cm"^) and shows that high-resolution observations yield higher deuteration values than single-dish 
observations. The NH2D observed transition, free of the optical depth problems that affect the NH3 lines in the core center, are a much 
better probe of the high-density nucleus and, in fact, its map peak at the dust continuum peak. Our analysis of the NH3 and NH2D 
kinematic fields shows a decrease of specific angular momentum from the large scales to the small scales. 

Key words. ISM: clouds - ISM: evolution - ISM: individual(L 1544) - ISM: molecules - Stars: formation - Techniques: interfero- 
metric 



1. Introduction est central density (IWard-Thompson et all 119991; iTafalla et all 



est central aensity QWara-inompsonetai.i ii^y^ iiaraiia et ai.l 
u u V. u ^ f , l2002HCrapsi et al.ll2005l) . It turns out that such cores have a few 

Dense starless cores are thought to be the forerunners of solar ^^^^^ properties in common such as high deuterium fractionation 
mass protostars and hence determining their physical and chem- ^ ^^^^^^ ^^^^ evidence for increased velocity disper- 

ical evolution is an essential step in our understanding of star ^-^^ indicative of contraction in the high density core nucleus, 
formation. Thus, measuring their physical characteristics (tem- Deuterium fractionation potentially gives us important informa- 
perature,density, velocity field) should give us insight into the ^^^^ ionization fraction at high densities while the ob- 

initial conditions just prior to the collapse of a protostar Most ^^^^^^ velocities and Une widths yield precious information on 
interesting in this regard are the cores of highest central density ^^g^^j^. momentum and infall. To exploit these to the full how- 
and column density which appear closest to the "pivotal struc- ^^^^ ^^^^^^^^ ^^^^^ ^^^^^^^ resolution measurements in a nearby 
ture" m which the core is on the verge of instability. ^^^^ ^^^^ ^^^^ i„ ^^^^ ^^^^^^^ ^^^^^^ ^^^^ q^^li^y 

Recent observations of both the dust continuum and molecu- terferometer maps of the nearby Taurus core L 1544 in the 23 
lar lines have allowed a selection to be made of the cores of high- GHz (J,K)=( 1 , 1 ) and (2,2) lines of ammonia (NH3)as well as in 

the 110 GHz 111 - loi transition of deuterated ammonia (para- 
NH2D). 

The choice of ammonia and NH2D for this study was partly 

Dascu oil ousci vauuiis caiiicu oui wiiii iiic iis.rtivi riaicau uc , . i i • i i r ^ ^ t^wrT~ri ^\ 

Bure Interferometer and the Very Large Array IRAM is supported by ^ased on the fact that it has been found (e.g. [Tafalla et al. | 

INSU/CNRS (France), MPG (Germany) and IGN (Spain). The National 2004) that nitrogen containing species such as ammonia and its 

Radio Astronomy Observatory is a facility of the National Science deuterium isotopologues trace preferentially the densest gas in 

Foundation operated under cooperative agreement by Associated nearby cores. The reason for this is not precisely known but is 

Universities, Inc. likely to be related to the relative volatility of both atomic and 
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molecular nitrogen. It is also important that the relative popu- 
lation of the (1,1) and (2,2) l evels of ammonia is an e xcellent 
temperature indicator (see e.g . lHo &~T ownes 1983; Danb y et"aLl 
ll988tlUngerechts et alJI 19861) due to the fact that radiative tran- 
sitions between the K-2 and K=l ladders of NH3 are forbid- 
den. We can therefore use the relative optical depths in these 
two lines as a measure of temperature. One of the basic ques- 
tions about pre-stellar cores is whethe r they h ave a positive tem- 
perature gradient outwards (see e.g. Leung 1975; Evans et aD 
l200ltlZucconi et aLllMTtrGani et al.l200'2[lLe'saffre et al.,.2005l) 



with temperatures of the order of 6 K in their center, as one 
would expect from theoretical considerations. On large scales 
> 0.05 pc, there is already evidenc e for higher dust temperatures 
on core edges than in core centers (fWard-Thompson et aLll2002l: 
iPagani et al.ll2004l: iKirk et al.ll2005T but on the scale of 0.01 pc 
(15" at 140 pc), the situation is less clear With this in mind, we 
obtained and present in this work D-array VLA NH3 observa- 
tions with a roughly 4"synthesised beam . We are thus sensitive 
to compact structure of size roughly 1000 AU. 

Measuring the temperature gradient is of importance mainly 
as a test of our understanding of the thermodynamics of pre- 
stellar cor es though it also i s of some importance for the core sta- 
bility (see lGalli et al.ll2002l) . The temperature measurements test 
our understanding of dust-gas coupling and CO cooling (see e.g. 
iBergin et al.l2006h . They are also important for the interpretation 
of dust continuum emission maps where it is usual to assume an 
isothermal dust temperature and spherical symmetry in order to 
derive the densit y distribution. The lat t er is often not v alid as 
we discuss later fconqalves et al.ll2004 iDotv et al.ll2005 b but is 
nevertheless a useful zero order approximation. Detailed mod- 
elling of L 1544 has been carried out by Doty et al. (2005) who 
fit the mm-submm continuum data of that source with a dust 
temperature drop from 15 K to 5.4 K. We note however that the 
ammonia single dish data for se veral sources including L 1544 
show no evidence for gradients (iTafalla et al.ll2002l 1200 4) and 
thus the existence of temperature gradients in pre-stellar cores is 
uncertain. 

L 1544 is a well-studied low-mass starless core in Taurus, 
which presents a large central density (> lO*" cm"^), consid- 
erable CO freeze-out (> 93%) and high degree of deuteration 
(N2D^/N2H^ ~0.25). These properties, combined with kinemat- 
ical probes of infall (differences in velocity between thin and 
thick tracers, and increasing line widths of thin tracers toward 
the core centre) makes this core a good candidate for being on 
the point of becoming unstable (Crap si et al.ll2005h . One expects 
the dynamics of the inner core to be better sampled by deuterated 
species and we have therefore obtained high resolution observa- 
tions of the In - loi transition of para-NHiDat 110 GHz with 
the IRAM Plateau de Bure interferometer (PdBI). The resolution 
is rather similar to that of the VLA and we are thus able to mea- 
sure the deuterium fractionation in the core nucleus. We also can 
derive important limits on angular momentum and infall. 

This paper is organized as follows, in section |2] we report 
the details of the VLA and PdBI observing runs, in section|3]we 
present our basic results, in section |4] we describe our analysis 
procedure and in section|5]we summarize our results. 

2. Observations 

We used the Very Large Array (VLA) to observe L 1544 in the 
(1,1) and (2,2) transitions of para-NHj (at 23694.4949 MHz and 
23722.6349 MHz respectively) and the Plateau de Bure inter- 
ferometer (PdBI) to observe the (lij - lo,i) transition of para- 
NH2D (at 1 10153.587 MHz). The reference position we used in 



thi s work is the dust continuum pe ak taken from the 1 .2 mm map 
of IWard-Thompsonetal.1 ^199% : (05:04:17.230, 25:10:47.70) 
J2000.0. 



2.1. VLA Observations 

Observations at the VLA were performed in compact configu- 
ration (D) on March 8"', 2003 and April 3''', 2003 for a total 
of 20 hours of observations. On-source observations were in- 
terleaved by the quasar 05 1 Oh- 180 (1 Jy at 1.3 cm according to 
the VLA calibrator manuaQ) for phase calibration, whereas ab- 
solute flux calibration was performed using 3C286 (2.59 Jy at 
1.3 cm). The difference in average flux density for 0510-1-180 
between the two observing runs was ~15% and we take this as 
the uncertainty in flux calibration. The NH3(1,1) and (2,2) data 
were acquired simultaneously using the 2 IF spectral line mode, 
sampling 64 channels over two 781 kHz (or 9.8 km s"') band- 
widths and reaching a channel width of 12.2 kHz (or 0.154 km 
s"'). While the (2,2) band was centred at the main hyperfine rest 
frequency, the (1,1) band was offset by 306 kHz to include both 
the main group of hyperfines and one secondary group (Fi =2- 
1 at 23693.905 MHz) of hyperfines. The VLA primary beam at 
the ammonia frequencies is 2', while the synthesized beam, after 
applying natural weighting, reached 4'.'34x3'.'45 with a position 
angle (PA) of 79^9. In order to increase the signal-to-noise ratio 
we reconstructed the clean map of the NH3(2,2) using a circular 
clean beam with full width at half power of 8". Calibration and 
data reduction were performed using the Astronomical Image 
Processing System (AIPS) of NRAO. 



2.2. IRAM Plateau de Bure Observations 

Plateau de Bure observations were performed in the C and D 
configurations (the most compact ones) on September 29''' 2002 
and December 18''' 2002 for a total of 16 hours of observing 
time. The correlator setup was chosen to simultaneously ob- 
serve the dust continuum at 110 GHZ and 230 GHz with a 
1 GHz bandwidth and the para-NH2D(l,l) and CO(2-l) with 
a resolution of 0.106 km s ' and 0.051 km s"' respectively. 
Phase calibration was performed using 0528-1-134 and 0415-1-379 
while absolute calibration was obtained using 3C454.3 for the 
D-configuration run and 3C273 for the C-configuration. The two 
phase calibrators have a 20% difference in flux in the two observ- 
ing runs and we take this as the calibration uncertainty. Cleaning 
was performed with uniform weight and the resulting synthe- 
sized beam is 5'.'8x4'.'5 in size and 8 1?4 in PA at 3 mm and 2'.'5 by 
2'.'2 with PA 81?0 at 1mm. The PdBI data were reduced and caH- 
brated using the Continuum and Line Interferometer Calibration 
(CLIC) developed by IRAM. 



3. Results 

3. 1 . Maps and spectra 
3.1.1. VLA Ammonia 

In figure [U we show the integrated intensity maps of NH3(1,1) 
and ( 2,2) from the VLA and the 100m antenna in Effelsberg 
(from lTafalla etai]|2002l FWHM ^ 37"). The qualitative differ- 
ence between the two is striking: while the intensity measured 
with the single dish follows very closely the column density 
traced by the dust emission, the VLA maps of NH3(1,1) peaks 
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20"south-east relative to the dust peak and the NH3(2,2) emis- 
sion shows a second peak also in the north west evocative of a 
"hole" in the emission. The angular size of the ammonia inte- 
grated emission seen with the VLA is 75"x36" in the (1,1) line 
and 68"x25" in the (2,2) line as compared to 105"x58" in the 
dust emission. 

In order to properly compare the interferometer and sin- 
gle dish maps, we cleaned the VLA data with a 40"beam to 
simulate the resolution of the 100-m data. The resulting map 
(Fig. [T]) shows that the two data sets are compatible. The dis- 
crepancy arises because the interferometer is insensitive to large 
scale emission and traces preferentially compact high density 
gas clumps. The peak intensities of these maps are a measure 
of the amount of flux lost in the VLA observations due to ex- 
tended structure: 68% in the NH3(1,1) (70% and 64% for the 
main hyperfine component and the first red hyperfine compo- 
nent, respectively) and more than 80% in the (2,2). As shown 
below the difference between these two is significant and thus 
the (2,2) emission appears to be more resolved out than the (1,1). 
This by itself is an indication that the temperature of the more ex- 
tended (presumably lower density) material in the L 1544 core 
is higher than in the compact nucleus. 

However, an important characteristic of our results is that the 
peak of the VLA (1,1) emission (effectively the ground state of 
para NH3) is clearly not coincident with the peak of the dust 
emission. We believe this is mainly due to the fact that the (1,1) 
emission is highly optically thick and thus our measurement to- 
wards the peak of the main (1,1) satellite reflects the temperature 
at the surface corresponding to unity optical depth. 

This is substantiated by Fig.|2]where we compare the spectra 
from the VLA with those observed with the single dishes. We 
measure at the dust peak a ratio of the satellite F=2-l emission 
to the main (F=l-1 and 2-2 blend) of 0.72 as compared to the 
expected ratio of 0.278 in the optically thin limit. This (see e.g. 
iHo & Towne"s ''l9"83) aflows an estimate of the main line optical 
depth and we derive a higher value from the VLA data (with 
an optical depth integrated over all the hyperfine components of 
T,oi > 30) than from the 100m (t,o, = 8.6), while 20" further 
from the peak the difference in opacity is reduced (tio, = 20 vs. 
Ttot = 7.7) and a smaller amount of flux is resolved out (we re- 
cover 75% of the flux). Our optical depths are in fact so high that 
they may affect the temperature gradients discussed later. These 
issues will be treated more quantitatively in section l43] where we 
report our radiative transfer calculations combined with the u, v 
modelling of the emission. For the moment, we note however 
that these results are consistent with a high ammonia column 
density in the neighbourhood of the dust peak. 

3.1.2. NH2D with Plateau de Bure 

Our NH2D data are extremely useful in that they confirm that 
ammonia remains in the gas phase at the densities of order 10'' 
cm"-' close to the dust peak. As in the case of ammonia, one can 
use the relative intensities of the hyperfine satellites to provide 
a measure of optical depth and in this case (see fig. |2| we find 
intensity ratios of the main component versus the first red com- 
ponent (=0.3) consistent with optically thin LTE (=0.278). Thus 
the fact that the NH2D intensity (see Fig. [T]i peaks at the dust 
peak can be considered a demonstration of the fact that NH2D 
is a tracer of the dense gas in the vicinity of the dust peak. It 
can thus be used to measure the kinematics in this region. The 
angular size of the NH2D emission seen by the interferometer is 
27"xl2" as compared to the dust emission size of 105"x58". 
Comparing the intensity at peak observed by the PdB with an 



IRAM-30m spectrum taken towards the nucleus of L 1544 we 
measure that 50% of the NH2D total flux is recovered by the 
interferometer. 

3.2. Kinematics 

We used the interferometric data in order to fit the NH3(1,1) and 
(2,2) and the NH2D spectra . In this manner, we obtained maps 
of velocity and line-width of the various tracers. Our results here 
are sensitive to differential flux loss in different channels as well 
as to the effects of opacity but, for the reasons discussed above, 
these problems are unlikely to be of great importance for NH2D. 
The central velocity information was used to comp ute the veloc- 
ity grad ient across the core, in analogous fashion to lCaselli et al.l 
(I20()2al) using N2H+ and N2D+. 

The velocity maps derived in this fashion are shown in Fig. [3] 
NH3(1,1) and (2,2) show the same general pattern with a to- 
tal gradient of ~9 km s ' pc"' with a PA of 160 degrees west 
of north. This pattern can be seen also in the channel maps 
presented in Fig [3] in both lines the emission is brighter in 
the north at blue velocities to become stronger in the south for 
higher velocities. Over a scale of 1 arc minute(0.04 parsec), this 
amounts to a specific angular momentum of 0.015 pc km s"' 
or 5 10^' c m^s~'. Our results m atch the velocity gradients cal- 
culated by ICaselli et al.l (l2002al) from single-dish observations 
N2D^, showing gradients of ~6 km s ' pc ' with a similar posi- 
tion angle as NH3. 

The different morphology of the NH2D maps compared to 
those of NH3 already tells us that the two species traces dif- 
ferent material, this is confirmed in the velocity maps, in fact, 
the NH2D velocity gradient map shows no clear sign of rotation 
with total gradient of ~2 km s ' pc ' with a PA of 60 degrees 
west of north. The coiTesponding specific angular momentum is 
0.001 pc km s"' or 3.5 10^" cm^s"', about an order of magnitude 
lower than the specific angular momentum measured with NH3. 
We note that the position angle obtained from NH2D is different 
from that seen in N2D^ by C aselli e t al. (2002a) and may suggest 
a difference between the bulk motions of neutral and ion species 
in the inner core. 

We note that the specific angular mo mentum obtaine d from 
NH3 is typical of starless cores (see e.g.. 'Ohashilll999l) . while 
the smaller specific angular momentum derived from NH2D is 
closer to the values of protostellar envelopes. Observations with 
higher velocity resolution should be needed to confirm this loss 
of specific angular momentum towards the small scales; how- 
ever, this differential rotational properties between the inner en- 
velopes and the ambient clou d is consistent with ot her observa- 
tions in more evolved objects jBelloche et al.ll2002l) . 

3.3. Dust Continuum and C0(2-1) 

The Plateau de Bure Interferometer observations, together with 
the NH2D(1,1) , simultaneously mapped L 1544 in the dust con- 
tinuum at 230 GHz and 110 GHz, and in CO(2-l). These ob- 
servations resulted in 3 o" non-detections of 0.9 mJy beam ' and 
0.5 mJy beam"' for the 1 mm and 3 mm continuum respectively. 
The corresponding limit for CO(2-l) was 80 mJy beam ' for the 
CO(2-l). 

Dust and CO em ission from L 1544 is easily detected by sin- 
gle dish telescopes dTafalla et al.| l2002'). but their distribution is 
clearly too extended even for the most compact PdBI configura- 
tions leading to a total loss of flux. The absence of compact con- 
tinuum emission at 1 mm confirms the starless nature of L 1544 
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and for an (arbitrarily) assumed temperature of 30 K yields an 
upper limit to the mass of a protoplanetary disk (on size scales of 
1 beam i.e. 2'.'5 by 2'.'2) of 1.8 ■ 10"'* solar masses. This estimate 
is a strong upper limit because we assume the small dust opac- 
ity for grains coated with thin i ce mantles and maximum size of 
0.25 yum) (equal to 0.78 cm^ g"' lOssenkopf & Henning|l994l) : if 
the maximum size of the dust grains in the disk grows to 1 mm 
the upper limit should be decreased by a factor of four (e.g. see 
dust opacities from Draine 2006). In fact, if L 1544 did contain 
a low-luminosity protostellar object si milar to that found b y the 
Spitzer Space Telescope in L 1014-IRS dYoung et al.l2004 . then 
we might have observed the continuum from the disk surround- 
ing the protostar (expecte d to be of few mJy in L 1014-IRS at 
1 mm. lYoung et al.ll2004l) . We assume here low optical depth 
conditions appropriate for large (> 100 AU) disks observed at 
millimeter wavelengths. For very small optically thick disks, our 
result amounts to a limit of Tdisk < 50K(Rdisk/lAU)"^. In any 
case we would have certainly observed a molecula r outflow sim- 
ilar to that in L 1014-IRS (which was detected bv lBourke et alj 
|2005| with a R.M.S. of 80 mJy (l"beam) ^ i.e. 4.8 times less 
sensitive than the present PdBI observations). 
Our data thus confirm the starless nature of L 1544. 

4. Analysis 

4.1. Evaluation of temperature in LIE 

Deriving temperatures from NH3 data has mostly been done 
using the so-called "standard analy sis" (e.g. lUngerechts et al.l 
Il986l: IStutzki & W innewisser' '1985V One can use for example 
equation A7 fromfUngerechts et al. ( 1986) to compute the "rota- 
tion temperature" T\2 defined by the ra tio of populations in then 
(2,2) and (1,1) levels. It is well known jWalmslev & UngerechtsI 
Il983h that at low temperatures this is in practice close to the ki- 
netic temperature (i.e. to LTE between the two levels). Thus, 

Tn = r2,2/ln[Tf//(0.45T^''2')] (D 

where 7^2,2 is the excitation energy of the (2,2) relative to the 
(1,1 level (41 K), and r'"' and t'^^ are the total optical depths of 
the (1,1) and (2,2) levels respectively (i.e the sum over all satel- 
lites). We assume here equality of line-widths and "amplitude 
factors" (see Ungerechts et al. 1986)) for the two lines. 

We thus see that our accuracy in determining the temperature 
depends on our accuracy in determining the optical depths. This 
is limited in the case of the ( 1 , 1 ) by the fact that using the relative 
intensities of the satellites becomes inaccurate for large optical 
depths (ratios close to unity). In the case of the (2,2) on the other 
hand, we do not measure the (2,2) satellites and thus are forced 
to derive the (2,2) optical depth from the ratio of the (2,2) and 

(1.1) main lines. Thus, we have 

T^2 = -I//22 In {1 - (Tm2/Tm,i) [1 - exp( -/ht'i"/)]} (2) 

where Tm,2 and Tm,i are the observed (2,2) and (1,1) brightness 
temperatures of the main hyperfine components and /22 and /i 1 
are the relative intensities of the main hyperfine structure for the 

(2.2) and (1,1) (0.5 and 0.8). The error-bars were evaluated prop- 
agating the errors on rM.2, ^m.! and Tj°'. Note that doubling the 
intensity ratio between the (2,2) and (1,1) transitions would re- 
sult in an increase of temperature of only 0.5 K, while doubling 
the opacity of the (1,1) transition would decrease the tempera- 
ture of 0.7 K. We expect that neither the errors on the opacity 
determinations nor the differential flux loss between the (1,1) 
and (2,2) reach this. Our main errors in this approach are caused 



by the breakdown of the isothermal assumption and to differ- 
ential flux loss between (2,2) and (1,1) (twice as much lost in 
(2,2)). 

In panel a of Fig.|4] we present the temperatures derived us- 
ing this technique as a function of projected distance from the 
dust peak. Both our pr esent VLA data set and the 100-m results 
of Tafalla et al.1 (|2002|) are shown. For comparison, we show re- 
sults from the theoretical model of Galli et al. (2002). A clear 
drop in temperature is found toward the continuum peak on the 
basis of the VLA data in contrast to the single dish data which 
show a rather flat profile. This is a clear suggestion that the com- 
pact high density core nucleus sampled with the interferometer 
is at a lower temperature than its surroundings. 

Our VLA data suggest that the gas temperature in the 
inner nucleus reaches values as low as 6 K. This contrasts 
with the minimum dust temperature of 11 K derived by 
IWard-Thompson et alj (!2()02) from ISOPHOT observations be- 
tween 90 and 200/im), although it should be noted that these ob- 
servations were done with ~80" of angular resolution and that 
far infrared observations alone are more sensitive to extended, 
warmer (T > 30 K) dust coming from smaller optical depths, 
like the outskirts of starless cores. In fact ou r observations are 
fully consistent with the best fit of iDotv et al.l (12005) who mod- 
elled the dust continuum at 450, 850 and 1300 jum) (hence typi- 
cal of cold dust) with a 3D radiative transfer code resulting in a 
central temperature of 5.4 K. 

In Fig. m we superimposed a modified version of t he the - 
oretical prediction (gre en curve) made b y iGalli et all ( |2002|) . 
Modifying the model in iGalli et all (|2002|) was necessary since 
this was calculated starting from a density profile that was de- 
rived from a constant dust- temperature of 8.75 K ( Tafal la et al.l 
I2OO2I) . Since the observed drop in gas temperature reflects a drop 
in dust temperature (the two are well coupled at hig h densities 
and low temperatures, see e.g. lGonqalves et all l2004l). we expect 
the central density to be enhanced using this new information. 
Since the temperature profile influences the estimate of the den- 
sity profile and vice-versa, we had to iterate few times the fol- 
lowing process; i) calculate the gas temperature from a given 
density profile ii) use that temperature profile to derive a new 
density profile fitting the 1.2 mm map iii) use this new density 
profile in point i). We converged after few iterations with the 
model represented in panel a) of Fig. |4] Note that this is not a 
fit of the temperat ure data but just an iterative process to make 
the model of Galli et al. I (I2002h self consistent with the 1.2 mm 
map. Nevertheless the agreement with the data is very good in 
the limits of validity of the model and the data. Including the 
temperature profile, the central density in L 1544 result to reach 
~ 2-10^ cm"^ while the size of the flattened area shrinks by 
30% and the slope of the density profile at large radii is not sig- 
nificantly changed. 

4.2. Column density and Deuteration in LTE 

Large deuterium enrichment is a characteristic of the high den- 
sity nuclei of pre-stellar cores and it is a useful test of our under- 
standing of the chemistry (e.g.'A ikawa et al.ll2005l: iFlower et al.l 
2006) to measure the [NH2D]/[NH3] ratio at high angular reso- 
lution. We have therefore combined our VLA and PdeB data for 
this purpose. 

We derived the column density of NH3 and NH2D in the 
constant excitation temperature approximation (CTEX, see also 
Caselli et al. 2002b). In the case of NH3, opacity and excita- 
tion temperature were derived through hyperfine structure fit- 
ting of the (1,1) transition using the formalism discussed earlier 
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We assumed an ortho-to-para ratio (o/p) of 1 when doing this 
and used a classical partition function at the temperatures of in- 
terest to correct for the excited states. The results from single- 
dish data can be co mpared to the column densities evaluated by 
iTafallaetalJ (12002 ) using the more sophisticated Monte Carlo 
approach. As shown in Fig 21 from the single dish data we de- 
rived a column density of NH3 of 3.5 ■ 10'"* c m~^, this , combined 
to a column density of H2 of 9.4 ■ 10^^ cm"^ dCrapsi et al. 2005). 
furnish es an abundance of 3.7 ■ 10 comparable to Tafalla et alj 
(|2002|) . The interferometric observations give a factor of ~2 
larger column densities, which may be an under-estimate in that 
our (1,1) optical depths are in some positions lower limits. 

In the case of NH2D, we are not observing the ground tran- 
sition - Oo,o at 332.8 GHz) and thus we need an estimate of 
the "rotation temperature". We solved this problem evaluating 

from th e ground state to th e obser ved level (1 1 1 - lo 1) using 
Eg. 2.21 of lWinnewisser et aP (Il978h : 

exp(-ro/r„) = expC-To/T'sB) X 

1 + (C,v/A,-,)exp(-ro/rt)[exp(ro/rfig) - i] 

l+(C,v/A,;)[l-exp(ro/rBB)] 

where Tq is the energy between the levels i and j, T^b is the 2.7 K 
background, the kinetic temperature, A,j is the Einstein coef- 
ficient and dj is the collision rate between the levels j and i. The 
latter can be expressed as C,j = cr,j (v,/,) n(H2), where cr,-^ is the 
cross section of the ; - i transition {- 8 ■ 10 cm"^, cf . iGreerJ 
1 19761; [Turner et al]| 19781) . (v,/,) is the thermal velocity dispersion 
and «(H2) the H2 volume density. Taking the kinetic temperature 
from the NH3 data and the «(H2) from the 1.2 mm data we were 
able to evaluate the excitation temperature of the ground tran- 
sition and thus the partition function of NH2D. This amounts 
to ignoring the splitting of the J=l level and is a crude proce- 
dure but we think it gives a reasonable estimate of the fraction 
of NH2D in the ground state. We assumed an o/p ratio of 3 for 
NH2D, consist ent with the average ob served NH2 D o/p ratio in 
starle ss cores jShah & Woottenll200lL but see also lFlower et alj 
I2OO6 I for a theoretical approach). 

The column density in the level {\\ \ - lo,i) was instead cal- 
culated in the CTEX limit. The profile of A^(NH2D) is shown in 
panel c) of Fig |4] In this case the column density derived from 
the interferometry data is only slightly higher than the estimate 
from the single dish. Dividing N(NH2D) by N(NH3) we were 
able to evaluate the degree of D-fractionation in NHsat the dust 
peak of L 1544 (see panel d) of Fig. |4|. The abundance ratio 
[NH2D1/[N H31 was foun d to be 0.5±0.2, thus higher than in 
N2H+ (0.26. ICasem etani20 02b). although caveats in the differ- 
ent flux loss between VLA and PdB apply. 

4.3. Model Fits of the Temperature and NH3 abundance 

In order to account for the temperature gradient along the line of 
sight as well as the complex flux loss from the interferometer, we 
made a model of the temperature, density and abundance struc- 
ture of the core. We then modelled the excitation of the molecule 
and the radiative transfer of the emission, and finally sampled the 
synthetic data cube in the u, v plane to compare with the VLA 
observations. 

The excitation and radiative transfer calculation were 
performed using the ID M onte Carlo radiative transfer code fo r 
molecular lines written by Hogerheiide & van derTakl(l2000l) . 
The code takes as input for each radial cell: the density of H2, 
the abundance of the molecule of interest, the dust and gas 
temperature and the radial velocity. Then, using the collisional 



rates from iDanbv et al.l ( Il988h and the energy levels, transition 
frequencies, Einstein A coefficients (from JPL j3 we can evaluate 
the population in each energy level and in each cell. This result 
is taken as input for the ray-tracing module to produce a data 
cube of the synthetic emission. We modified the ray-tracing 
module to include the hyperfine structure of NH3. The output 
data cube can then be sampled in the u, v plane to reproduce 
the VLA observations. We produced such synthetic visibility 
datasets with the MIRIAD task UVMODEL. The Monte Carlo 
co de used in thi s paper was carefully tested against the one used 
bv lTafalla et al.l (|2002|) resulting in a very good agreement. 

We parametrized the structures of density, temperature and 
abundance in the following way: 



«(r) = 



no 



1 + (r/ro,,,)"" 



T(r) = To, ~ 



1 + ir/rojy 



The forms for the density and abun dance profiles where cho- 
sen for consistency with the results of iTafalla et alj (l2002l ). The 
form chosen for the temperature structure has main ly the virtue 
of sim plicity. It is consistent with those computed bv lEvans et al.l 
(I2OOII) (e.g. that shown in their Fig. 6 can be represented by 
Tout = 13.1 K, Tin - 6.5 K and roj = 26") though we note 
that such models neglect the breakdown of gas-grain coupling in 
the outer parts of the core. We therefore ran the calculation for 
a grid of values for Tom, Ti„, ro.r, Xq, ax and compared the re- 
sults with the observations of the NH3(1,1) main and first hyper- 
fine components and the NH3(2,2) main component observed at 
VLA and Effelsberg. The density distribution is not constrained 
by the NH3 observations but by the 1 .3 mm dust continuum. For 
each temperature distribution we then recomputed the density 
distribution that would fit the 1.3 mm observations in similar 
fashion to Tafafla et al. (2002). 

The quantities we wish to match with both the VLA and the 
Eff'elberg observations are 

- the integrated amplitude vs. u, v distance profile for each of 
the components (NH3(1,1) main and first hyperfine (hf) and 
NH3(2,2) main) 

- the profile of the ratio between NH3(2,2)/ NH3(1,1)/,/ 
(which roughly speaking is a measure of the temperature) 
and the ratio between NH3(l,l),„a,„/ NH3(1,1)/,/ (which 
roughly speaking is a measure of the column density) 

- the line profile of all the observed data 

We compared the amplitude profile both versus the u, v distance 
of the baselines and versus the distance to the centre of the 
source in the inverted imag^ Those plots are equivalent since 
they differ only for a Fourier transform, however since the plot 
of amplitude versus distance from the source centre shows bet- 
ter the emission in the centre of the source we choose to present 
only those plots. 

The model that was found to best fit all this quantities has 
Tout = 12 K, Ti„ = 5.5 K, ro.r = 18", = 8 ■ lQ-\ ax = 0.16. 



- these informations are collected in the Leiden Atomic and 
Molecular Database at http://www.strw.leidenuniv.nl/~moldata/ 

' for this comparison we did not "clean" either the data or the model 
since this would have introduced further uncertainties; comparing the 
"dirty" maps is instead exactly equivalent to the comparison in the u, v- 
amplitude plane. 
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The density, temperature and abundance profiles of this model 
are shown in fig.|5]in solid lines. As a comparison we show with 
dashed lines the profiles relative to the best fit of the single -dish 
data alone with constant temperature dTafalla et alj|2002l) . The 
new temperature profile is warmer in the outside and colder in 
the very centre and this results in a more peaked density struc- 
ture (no -2.1 10^ cm"-', ro „ = 14", an = 2.5). Th e abundance 
rise in the centre of the core shown bv lTafallaetan (.2002) IS re- 
quired also in our best fit although it is less steep (ax = 0.16 vs 
0.3). 

The reason for these changes in the temperature structure can be 
seen in the comparison between the observations and the pre- 
dictions of the two models presented in Fig |6] In each panel of 
Fig|6]the data are shown in grey, the black-bold curves represent 
the prediction of our best fit model while the thin-red lines show 
the best fit of the single dish data assuming constant temperature. 
The lower temperature in the centre is required to fit the intensity 
of the NH3(2,2) line which is highly over produced in the con- 
stant temperature model (see panels g and 1). Lowering the cen- 
tral temperature requires the external temperature to be increased 
in order to still fit the single dish observations of NH3(2,2) (panel 
c). The quality of this new temperature profile can be checked 
also against the ratio of the NH3(2,2) and the NH3(1,1)/,/ which 
is very sensitive to variations in the temperature profile (panel h). 
The constant temperature predictions are dramatically different 
from the observed ratio of NH3(2,2)/NH3(1,1);,/. Adjustments to 
the abundance profile are finally needed to fit the intensity pro- 
file of the two observed NH3(1,1) hyperfine components. The 
characteristic profile of the NH3( 1,1 ),„«,„ (panel e) and the flux 
loss of the central channels of the main component (panel k) 
are two sensitive indicators of a good fit to the abundance pro- 
file. In panel d the ratio between the two NH3(1,1) hyperfine 
components shows the overall quaht y of the adopt ed abundance 
profile. We note that the best fit of lTafallaet all (^002) better 
reproduces the intensity profile of the first red hyperfine com- 
ponent of NH3(1,1) observed with the single-dish (panel b and 
i). However increasing the intensity of the NH3(l,l)fty would re- 
quire higher abundances or warmer inner temperatures and this 
would significantly worsen the modelling of the interferometric 
data (see panels d, f and h). A possible solution of the problem 
could be found introducing a 2D geome try (see e.g. the rota- 
tionally flattened envelopes of lUlrichlll976 i) but, since the effect 
on the temperature profile would be modest, this is beyond the 
scope of the present paper. 

5. Conclusions 

We presented high-resolution, interferometric observations of 
NH3 and NH2D toward L 1544, a low-mass starless core on the 
verge of forming a star The gas temperature and deuterium frac- 
tionation are derived for the first time with a resolution smaller 
than 800 AU. The main results of this work are: 
1) The gas temperature drops to ~5.5 K at the cor e density peak, 
in agreement with theoretical model pred ictions (IZucconi et al.l 



l200lt lEvans et al.l 



uum modelling of 



2OOU iGal i et al.ll2002h and wifli flie contin- 



Dotv et al 



(|2005). In particular, our data can 
be reproduced by the models of Galli et al. (2002), simultane- 
ously fitting the dust continuum data. A similar con clusion about 
the ga s temperature drop has been obtained by iPagani et al.l 
(120071) based on N2H+, N2D+, and NH3 observations of L 183. 

2) The central density derived from literature data but including 
the temperature profile reaches 2 ■ 10^ cm"^, 50% higher than 
considering constant temperature at 8.75 K. 

3) The column densities for both NH3 and NH2D are (factor 



of ~2) larger than previously found with s ingle-dish antennas 
(iTafa Ua et al. l l2002t IShah & WoottenI 1200 l l). evidence that we 
are probing a higher density gas and that, unlike CO, the am- 
monia freeze-out is not significant in the core centre, where 
niHjJ ^ 10^ cm-\ 

4) The deuterium fractionation in NH3 is 0.5+0.2, almost con- 
stant in the central 4000 AU, where NH2D is detected 

5) The NH2D map, unlike those of NH3, peaks at the dust con- 
tinuum peak showing that it is a very good tracer of the high den- 
sity gas. A comparison between the velocity field extracted from 
NH2D with the one derived from NH3 shows an order of mag- 
nitude decrease of specific angular momentum from the large 
scales to the small scales. 
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40" res. 



Dust peak position {0",0"] 



NH3(2,2) peak position (10", -20"; 




- 6200 AU 




Ka [J3000] 

Fig. 1. Integrated intensity maps of L 1544 in NH3(1,1) and 
(2,2) and para-NH2D(ln - Iqi). Dashed contours in all the 
panels are the 50% and 95% levels of the 1.2nim continuum 
taken from I Ward-Thompson et al] (1 19991) and smoothed to 22". 
Left column: Observations taken with VLA or PdBI were in- 
tegrated over all the hyperfines available (NH3(1,1) values re- 
fer only to the main component and the first red component). 
Angular resolution is 4'.'34x3'.'45, 8"and 5'.'8x4'.'5 for NH3(1,1), 
NH3(2,2) and NH2D(1-1), respectively, and it is reported in the 
bottom right. Levels are 15% to 95% by 20% of the map peak 
which is 0.25, 0.05 and 7.0 K km s"' for NH3(1,1), NH3(2,2) 
and NH2D(1-1), respectively. Right column: Integrated inten- 
sity maps of the observations from the single dish (top) and in- 
terferometer data smoothed to 40" (bottom). Levels are spaced 
percentually as in the left column but peak values are 8.5 K km 
s 'and 2.8 K km s 'for the single dish map of NH3(1,1), and the 
smoothed VLA map of NH3(1,1), and 0.47 Kkm s 'and 0.08 K 
km s^'for the NH3(2,2) maps from Effelsberg and VLA. A 40" 
beam is reported at the top-left corner. 
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Fig. 2. Comparison between NH3 ( 1 , 1 ) and (2,2) and NH2D( 1 - 1 ) 
spectra observed with the interferometers (VLA and PdBI) and 
the same line observed with the single-dish (100m in Effelsberg 
and IRAM-30m). The three left box show the spectra observed 
toward the dust emission peak while the two right ones were 
taken toward the NH3(2,2) peak (offset by 10",-20"). 
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Fig. 5. Temperature, density and NH3 abundance profile that best 
fit the interferometric and the single dish observations simultane- 
ously (solid lines). The profiles relative to the best fit to the sin- 
gle dish data alone bv Tafalla et aL (,200Z) are reported in dashed 
lines for comparison. 
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Fig. 3. Kinematics in L 1544. Top row: velocity gradients calculated in 9 adjacent points are overlaid on the integrated intensity 
map for each observed transition. The direction of the arrow points to the increasing velocity while the lenght measure the magnitude 
of the gradient (a 10" arrow corresponds to a gradient of 33 km s ' pc The thick arrows in the bottom-right represent the total 
gradient evaluated with all the fits (the magnitude for those is ten times smaller: i.e. a 10" arrow corresponds to a gradient of 3.3 km 
s"^ pc"^) Bottom row: channel maps of the main component of each observed Une (VLA and PdBI). 
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Fig. 4. Temperature and column densities profiles as derived in the LTE approximation (see sections 143] and Wl% . Observations 
from VLA are represented wi th empty squares, while those obtained at the 100m single-dish in Effelsberg are shown with filled 
triangles. Literature data from 'S hah & WoottenI (1200 lb are shown with filled circles. Panel a): Gas temperatur e profile in L 1544 
both from single-dish (triangles and circle) and VLA data (empty squares). The solid line is an adaptation of the lGalli et alj ( |2002[) 
model now using a central density of 2 ■ 10^ cm""* which is consistent with the low temperature observed (see section [431 ). Panel 
b): NH3 column density profile from both single dish data (triangles and circle) and VLA data (empty squares). Panel c): same as 
panel b) but for NH2D. Panel d): A^[NH2D]/A^[NH3] profile for the literature data and the new high-resolution data. 
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Fig. 6. Comparison between the observations of NH3 (1,1) and (2,2) taken with the 1 00m antenna of Effelsberg and the VLA, and the 
model prediction. Data are always displayed in light grey. The best fit is shown in bold black lines, while the thin red curves show 
the best fit of the single dish data alone. In panel a) we show the integrated intensity of the main hyperfine component of NH3(1,1) 
as observed with the single-dish. Panel b) and c) are similar to panel a) but for the first red hyperfine component of NH3( 1 , 1 ) and the 
main component of the NH3(2,2). Panels d) to h) present the intensity profiles observed with the VLA compared with the models. 
Both the data and the models were inverted from the u, v plane to the image plane using a fast fourier transform and then azimuthally 
averaged; the consequent negative values are not due to errors but to the negative side-lobes, and we chose to include them in the 
modelling rather than add further uncertainties due to cleaning. In panel d), we show the intensity ratio profile between the first 
red hyperfine component of NH3(1,1) and the main component, this is an indicator of column density distribution; the asymptotic 
behaviour at 25" and 60" is due to the diff'erent extension of the two maps that determine a different size of the negative side-lobes. 
In panel e), f) and g) we present the intensity profiles of the NH3(1,1) main component, NH3(1,1) first red hyperfine and NH3(2,2) 
main component, respectively. Panel h) shows the amplitude ratio profile between the NH3(2,2) main component and the first red 
hyperfine component of NH3(1,1), this is an indicator of temperature. From panels i) to 1) we compare the spectra at peak observed 
with the single dish and the interferometer with the models. 



